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ABSTRACT 

It is now well established that Kepler's supernova remnant is the result of a Type la explosion. With 
an age of 407 years, and an angular diameter of ~ 4', Kepler is estimated to be between 3.0 and 7.0 
kpc distant. Unlike other Galactic Type la supernova remnants such as Tycho and SN 1006, and SNR 
0509-67.5 in the Large Magellanic Cloud, Kepler shows evidence for a strong circumstellar interaction. 
A bowshock structure in the north is thought to originate from the motion of a mass-losing system 
through the interstellar medium prior to the supernova. We present results of hydrodynamical and 
spectral modeling aimed at constraining the circumstellar environment of the system and the amount 
of 56 Ni produced in the explosion. Using models that contain either 0.3M Q (subenergetic) or 1.0M Q 
(energetic) of 56 Ni, we simulate the interaction between supernova la ejecta and various circumstellar 
density models. Based on dynamical considerations alone, we find that the subenergetic models favor 
a distance to the SNR of < 6.4 kpc, while the model that produces 1M Q of 56 Ni requires a distance 
to the SNR of > 7 kpc. The X-ray spectrum is consistent with an explosion that produced ~ 1M Q 
of 56 Ni, ruling out the subenergetic models, and suggesting that Kepler's SNR was a SN 1991T-like 
event. Additionally, the X-ray spectrum rules out a pure r~ 2 wind profile expected from isotropic 
mass loss up to the time of the supernova. Introducing a small cavity around the progenitor system 
results in modeled X-ray spectra that are consistent with the observed spectrum. If a wind shaped 
circumstellar environment is necessary to explain the dynamics and X-ray emission from the shocked 
ejecta in Kepler's SNR, then we require that the distance to the remnant be greater than 7 kpc. 
Subject headings: Hydrodynamics, ISM: individual (SN 1604), Nuclear Reactions, Nucleosynthesis, 

Abundances, ISM: Supernova Remnants, Stars: Supernovae: General, X-Rays: 

ISM 



1. INTRODUCTION 

Type la supernovae (SNe) are believed to be the 
thermonuclear explosion of a C+O white dwarf (WD) 
that is destabilized when its mass approaches the Chan- 
drasekhar limit by accretion of material from a binary 
companion. After the central regions ignite, a burn- 
ing front propagates outwards, consuming the entire star 
and leading to a characteristic ejecta structure with ~ 
0.7 M of Fe-peak nuclei (primarily in the form of 56 Ni 
that powers the SN light curve) in the inside, and an 
equal amount of intermediate mass elements on the out- 
side. While the amount of 56 Ni can vary between ~ 
0.3 M© and ~ 1 M Q , the stratification of the ejecta re- 
main s constant across SNe la subtypes (Mazzali et al. 
[2007[). The simple structure of the ejecta results in uni- 
form light curves and spectra and make SN la useful dis- 
tance indicators (e.g. JRiess et al.lH998tlPerlmutter et al.l 
1999). Depending upon the nature of the WD compan- 
ion, SN l a progenitors are di vided into either single de- 
generate (IHachisu et al J 1996ft or doub le degenerate cases 
(|Iben fc TutukovI1984 Wcb binkil 19841) . Finding clear ev- 
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idence supporting one of these two scenarios has become 
a cornerstone in stellar astrophysics, as unknown evolu- 
tionary effects associated with the progenitors may intro- 
duce systematic trends that could limit th e precision of 
cosmological measurements based on SN la ([Howell et al.1 
120011) . 

Most Type la supernova remnants (SNRs) reflect 
the simplicity of the explosion. Galactic SNRs 
such as T ycho and SN 100 6 appear remarkably 
symm etric ([Warren et al.l 120051 : ICassam-C henai' et al.1 
2008), unlik e their core-collapse counterparts , such 
as Cas A (iDeLanev et al.l 120101: lLopez et al.l 120111: 
iHwang fc Lamind I2012D or younger extragalactic SNe 
such as SN 1980K or SN 1993J ([Milisavlievic et al.1 
2012). The morphology and X-ray spectra of SNRs 
with known ages and a Type la classification also indi- 
cate that the progenitors do not modify their surround- 
ings in a strong way - in particular, while absorption 
features seen in the spectra of Typ e la SNe suggest 
small cavities of ra dius > 10 17 cm dPatat et alJ l2007t 
Blondin etaD 120091: iBorkowski et~aTl l2009t iSimon et all 
20091 iSternberg et al.l 120111 ) , there is little evidence for 
large (r ~ 3 - 30 pc) wind-blown ca vities around 
the explosion site (Badenes et ail I2007D . However, 



iWilliams et al.l ([201 ID recently showed that the candi- 
date Type la RCW 86 is currently interacting with a cav- 
ity wall and that the cavity has a radius of ~ 12 pc, in the 
middle of the range of sizes expect ed from the accretion 
wind models ([Badenes et al.ll2007[ ). The lack of a size- 
able population of la SNR inside of cavities is particularly 
relevant as the current paradigm for SN la progenitors 
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in the single degenerate channel predicts fast, optically 
thick outflows from the W D surface which wo uld leave 
behind such cavities (e.g. lHachisu et al.l Il996j ). Addi- 
tionally, the dynamical and spectral properties of young 
la SNRs like Tycho and 0509-67.5 are more consistent 
with an interac tion with a constant densi ty interstellar 
medium (ISM) (IBadenes et al.l[2006l [2008h . 

A notable exception to this is Kepler's SNR (G4.5+6.8; 
Kepler), shown in X-rays in Figure [T] Kepler's SNR 
has recently been firmly established as the remnant of 
a Type la SN, b ased on the O/Fe rat io observed in the 
X-ray spectrum (jRevnolds et al.l [20071 ). However, unlike 
many other Type la SNRs, the morphology of this « 400 
yr old remnant shows a strong bilateral symmetry, and 
optical observations reveal clear signs of an interaction 
between the blastwave and a dense, nitrog en-enhance d 
circumstellar shell with a mass of - IMq (IVml[2008h . 
Near infrared observations also show evidence for circum- 
stellar kno ts akin to the quasi-stat ionary floculi of Cas- 
siopeia A (|Gerardv fe Fesenl I2001D . This suggests that 
either the progenitor of Kepler's SN or its binary com- 
panion might have been relatively massive, creating a 
bowshock shaped circumstellar medium (CSM) structure 
as the system lost mass and moved against the surround- 
ing ISM. lBorkowski et al.l (|1992D found values of M ~ 5x 
10~ 6 M Q yr -1 and v W i n d = 15 km s _1 are required in 
order to explain the morphology of Kepler's bowshock. 
While it is likely that not all progenitor scenarios for 
SNe la are the same, the dense CSM wind required to 
create the observed morphology appears to be in conflict 
with o bservations of the r e cent nearby SN la 201 lfe in 
M101 ([Nugent et al.|[20TTI: iLi et al.ll2CKLl"l) . In the case 
of 201 lfe, the radio and X-ray observations place con- 
straints on the circumstellar environment with ncsM ^ 
6 cm -3 , or in the case of a st ellar wind, M < 6xl 0~ 10 
M yr" 1 (jHoresh et al.l [20121 : [Chomiuk et aljl20"H . It 
is worth noting, however, that the recent observations 
of SN 2011fe probe spatial scales of ~ 10 16 cm, corre- 
sponding to the last ~ 100-1000 yr of the progenitor's 
evolution, while the blastwave of Kepler's SNR is probing 
the mass-loss history 10 5 - 10 6 yr ago. 

Located well above the Galactic plane, the SNR has 
an angular s ize of 2' in radius, but t he distance is poorly 
constrained. iRevnoso fe Gossl (|1999) used H I absorption 
to estimate a distan ce of 4.8 < Dsnr < 6.4 kpc, while 
ISankrit et al.l (|2005f) estimate that the SNR is somewhat 
closer, with a lower limit of 3.0 kpc on the distance, based 
on the width of the Ha line and proper motion mea- 
surements. The non-detection of TeV 7-rays by HESS 
(F 7 < 8.6 x 10 -13 erg cm~ 2 s _1 ) pla ces a lower limit 
on th e distance to Kepler of ~ 6 kpc (jAharonian et al.l 
[20081) assuming a mo del for the gamma-ray emission 
([Berezhko et al.l 120061 ). for a normal la SN. To match 
the angular size, the distance would have to be even 
greater than 6 kpc (> 7 kpc) if the SN were more ener- 
getic than a normal la, or else the blastwave would have 
to exp and into a CSM with considerable density. iVinkl 
( 2008) measured the velocity of the SNR forward shock 
and estimated v s = 4200(i4kpc km s _1 and concluded that 
the distance must be > 6 kpc or otherwise the SN was 
subenergetic. The distance estimates are incompatible 
with one another: either 3.0 < Dsnr < 6.4 kpc, or Dsnr 
> 7 kpc, based on the non-detection of TeV gamma-rays, 



an important new constraint. Given the uncertainty on 
the distance, the blastwave has a radius between 1.75 > 
Rfs > 3.7 pc, or Rps > 4.1 pc. 

Here we present hydrodynamical simulations of Type 
la SN ejecta interacting with an external medium shaped 
by a circumstellar wind. Combining models for the pre- 
supernova mass-loss with models for the ejecta density 
distribution and composition, we simultaneously fit both 
the morphology of the SNR and the observed X-ray spec- 
trum within the observable constraints- the age and an- 
gular size of the SNR along with the bulk spectral proper- 
ties. Since the circumstellar environment is complicated 
by the bowshock structure in the north, we focus our at- 
tention on the southern portion of the SNR (see Fig. [IJ , 
in an attempt to minimize the number of free parameters 
in our hydrodynamical simulations. Based on our joint 
hydrodynamical and spectral fits, we conclude that Ke- 
pler was a luminous Type la SN. If the SNR is expanding 
into a stellar wind we require a distance to Kepler of > 
7 kpc. If, on the other hand, the SNR is expanding into 
a uniform medium, the distance is somewhat closer, ~ 5 
kpc. 

2. X-RAY EMISSION FROM KEPLER'S SNR 

A 750ks Chandra image of Kepler's SNR (PI: 
S. Reynolds) is shown in Figure [TJ with the X-ray emis- 
sion from the pie-sliced region shown in Figure [2] (left). 
As shown in Fig [2] (left), the X-ray spectrum is domi- 
nated by emission from Fe K and Fe L, as well as emission 
from intermediate mass elements such as silicon, sulfur, 
and argon. The amount of iron emission in the spec- 
trum is a key constraint on the explosion, as it will al- 
low us to discriminat e between Type la explosion models 
(jBadenes et al.l[2~005l ). Besides emission from the shocked 
ejecta, there is also emission from circumstellar material- 
mainly oxygen and neon, which over laps with the L-shell 
iron emission (jRevnolds et alJfeOOTl ). as well as nonther- 
mal continu um which contr i butes to the emission above 
a few keV (jRevnolds et al.ll2007t iCassam-Chena'i et al.l 
l2008f l. In the subsequent discussion on the fits to the 
spectrum, we do not account for either the contribution 
from the swept up circumstellar material or nonther- 
mal emission in the spectral models- the contribution 
by shocked circumstellar material to the thermal X-ray 
emission is small a nd accounts for only a few percent of 
the total emission ([Reynolds et a l. 2007), while the non- 
thermal emission serves to raise the continuum above ~ 4 
keV. The processes that lead to th e nonthermal emission 
could alter the thermal spectrum (jPatnaude et al.ll2009l 
I2010D . but those affects h ave been shown to be small, par- 
ticularly in young SNRs (jPatnaude et alj 2010). Interest- 
ing characteristics of the spectrum include: the presence 
or absence of emission from iron, which will constrain 
the explosion energetics; the centroids of K-shell emis- 
sion lines, which constrains the ionization state of the 
shocked gas; and flux ratios of individual lines, such as 
Ka to H-like Lya, which constrain the density of the 
circumstellar material. 

The spectra from the combined 750ks observation 
shown in Figured] was extracted from the Level 2 Chan- 
dra event list using CIAO 4.1, and weighted response 
matrices were computed using CalDB version 4.1.3. The 
broadband spectrum (Fig. [2] left) shows significant emis - 
sion from both Fe L and Fe K ([Reynolds et al.l 12007). 
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Given the complexity of the Fe L emission at ~ 1 keV, we 
focus only on those emission lines above 1.5 keV, as they 
are well separated from neighboring lines. To measure 
the emission line centroids, we modeled the spectrum 
as a powerlaw continuum with a series of Gaussians for 
each emission line. The spectral fit is shown in Figure [5] 
(right) . We chose to model t he con tinuum as a powerlaw 
since iCassam-Chena'i et al.l f|2008T ) showed that much of 
the continuum emission in Kepler is nonthermal in na - 
ture (see Figs. 6 and 8 of lCassa~ Chenaiet alJl2008h . 
For the spectral fits, we used XSPEC Version 11.3.2ag. 
We used a region nearby the supernova remnant for a 
local background subtraction. Since we are only fitting 
the data above ~ 1. 5 keV, we do not inclu de absorption 
in our spectral fits. iRevnolds et al.l (|2007l) found an ab- 
sorbing column of N# = 5.2 x 10 21 cm -2 by fitting fea- 
tureless regions along the remnant's rim to an absorbed 
srcut model. Fitting only the data above 1.5 keV, we 
were unable to constrain the column density. However 
in the subsequent sections where we compare our hydro- 
dynamical and spectral models to the data, we assume a 
column of 5.2 x 10 21 cm~ 2 so that we may compare the 
computed X-ray emission below 1.5 keV directly to the 
data. 

The results of our centroid fits are listed in Table [TJ 
The key results in Table [1] are the line centroids, and 
in the subsequent discussion, we will focus on the line 
centroids of of Si Kq, S Kq, and Fe Ka. We also point 
out in Figure [2] (right) that there is little emission at 2 
keV from Si XIV Lya. This is confirmed in our spectral 
analysis where we are only able to set an upper limit 
on the Si XIV Lya emission from shocked ejecta in the 
south. The lack of emission at this energy will be used 
to constrain the CSM structure and density. 

3. HYDRO DYNAMICAL MODELING 

In this section, we present results from a modeling 
effort where we have coupled the hydrodynamics to a 
nonequilibrium ionization (NEI) calculation to produce 
joint spectral and dynamical fits to Kepler's SNR. We 
describe, in general terms in § 3.1 the formation of the 
bowshock structure in the north of Kepler and how it can 
be used to constrain the mass-loss parameters of the pro- 
genitor system. We also describe how the choice of la ex- 
plosion model (e.g. the explosion energy and amount of 
56 Ni that is synthesized) imposes a constraint on the un- 
certain distance to the SNR. In the subsequent sections, 
we compare the bulk spectral properties of Kepler's SNR 
to those predicted by our spectral models. 

3.1. Morphology of Kepler's SNR in the Context of a 
CSM Interaction 

The presence of a bowshock and shell of circumstel- 
lar material in the north of Kepler's SNR suggests that 
there was an outflow from the progenitor or its compan- 
ion prior to the explosion. This is because the bowshock 
can only form at a converging flow, in this case a wind 
interacting with the interstellar medium, where the ISM 
has a net inflow velocity relative to the progenitor, due 
to the progenitor's sy stemic velocity through the ISM. 
iChiotellis et al.l (|2012l ) recently modeled the dynamics of 
the SNR blastwave through the bowshock and concluded 
that the dynamics of the SNR blastwave are consistent 



with the presence of a bowshock. We can use estimates 
for the location of the bowshock stagnation point as a 
constraint on the wind parameters. Generally, the lo- 
cation of the stagnation point, r n , is w ritten as (e.g. 
iHuang fe Weigertill982UBorkowski et al.lll992fl : 

/ • \ 1/2 

1 to ir>3 / Mv wind \ ,,, 

r = 1.78 x KT 7 pc , (1) 

where M is the mass-loss rate, i\vmd is the wind speed, 
!ia, n b is the ambient ISM density, and is the systemic 
velocity of the star. In order to explain the large dis- 
tance above the Galactic plane, the progenitor likely had 
a lar ge systemic velocity, u+ ~ 250 km s _1 (jBandieral 
1987), a notion confirmed by the high proper motion and 
radial velocities of nitrogen-rich knots in the remnant, as 
well as the narro w component of the Ha line in the non- 
radiative shock. iBorkowski et al.l (|1992l ) estimated that 
the undisturbed ISM material had a density of 10~ 3 - 
10 -4 cm" 3 , corresponding to the hot component of the 
ISM (|McKee fe OstrikeHll977l ). Equation Q] can be in- 
verted to estimate the mass-loss rate of the wind: 

M = 3.15 x w -r?»E^ (n\ . ( 2 ) 

"wind \P C / 

Given the uncertainty of greater than a factor of two on 
the distance, the stagnation point is loca ted at a radius of 
2-4 pc , similar to the 2-3 pc assumed by IBorkowsk i et all 
(1 19921). The amount of mass in the shell is ~ 1 M©. 
IChiotellis et al.l (|2012|) have suggested that the nitrogen- 
rich knots in the CSM, combined with the solar mass of 
material in the shell, argue in favor of the CSM being 
sculpted by the outflow from an asymptotic giant branch 
(AGB) star. The mass-loss rates from AGB stars can 
vary over several orders of magnitude , but the velocity of 
the w ind is generally 10-2 km s" 1 (pchoier fe Olofssonl 
f200ll lOlofsson et al.il200l iRamstedt et al.ll2009| ). Using 
this range in velocity combined with the range in the 
radius of the stagnation point yields mass-loss rates of 
10~ 6 - 10" 5 Mq yr- 1 . 

To model a CSM wind, the mass-loss rate and wind ve- 
locity are all that are required, since p a mb = Ar~ 2 where 
the wind normalization A — M /(An u w i n d)- We couple 
this model for the circumste llar environment to m odels 
for Type la SN ejecta from iBadenes et al.1 ([2003). We 
choose to focus on models DDTa (Esn = 1.4xl0 51 erg) 
and DDTg (E SN = 0.9xl0 51 erg). The choice of models 
is designed to encom pass the diversity of la explosions 
(jBadenes et al.l 120081 ) . In a constant density environ- 
ment, the choice of explosion model does not significantly 
influence the radius and velocity of the forward shock as 
long as the kinetic energy is ~ 10 51 erg. In a wind envi- 
ronment, the difference in forward shock radius between 
an energetic and subenergetic model, where the ejecta is 
approximated as a powerlaw in velocity, is a little larger, 

since R FS cx Eg(j ( " 2) rather than Eg^ n (|Chevalierlll9jl . 
Of more interest is the resultant X-ray emission, which 
is influenced by the density and chemical composition of 
the ejecta. For instance, in the DDTa model, ~ 1M© 
of 56 Ni is produced, compared to 0.3M© in the DDTg 
model. Since the 56 Ni decay is responsible for the iron 
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observed in the ejecta, the amount of 56 Ni produced in 
the SN will profoundly affect the emitted X-ray spec- 
trum. 

Using the VH-1 hydrodynamics code, a numerical hy- 
drodynamics code developed at the University of Vir- 
ginia by J. Hawley, J. B londin, and collaborators (e.g., 
iBlondin fe Lufkinl [19! )3), we modeled the expansion of 
the SN ejecta into the CSM wind to an age of 400 yr. 
We varied the wind speed between 10-20 km s" 1 and 
the mass-loss rate as described above. The parameter 
space results in a grid of models that relate the radius of 
the forward shock with the wind mechanical luminosity, 
L wind = ^M^ ind , for each explosion model. In Figure [3] 
we plot the shock radius as a function of wind luminosity 
for each combination of wind speed and mass-loss rate. 

Roughly speaking R FS oc L^ d cx v^ d , for a given mass- 
loss rate. This is because the CSM density normalization 
is set by the wind speed - higher wind velocities transport 
the mass further away from the progenitor, while lower 
velocity winds result in more mass deposited closer to 
the progenitor. In the higher wind velocity models, the 
blastwave propagates further after 400 years since it has 
to move through less material. 

As discussed earlier, the distance to Kepler is not well 
known, ranging from 3 - 6.4 kpc, or > 7 kpc. In FigureE] 
we also note the distance to Kepler on the right hand 
axis, using the measured angular ra dius of 2'. The upper 
limit of 6.4 kpc ([Sankrit et al.ll2005j ) on the distance rules 
out all DDTa models on dynamical considerations alone. 
However, if the distance is > 7 kpc, as expected from the 
TeV gamma-ray non-detection, then most DDTa models 
are allowed, while a sizeable fraction of DDTg models 
(M > 6x 10~ 6 M Q yr -1 ) are no longer favored. 

3.2. Spectral Fitting: Wind Models 

To compare our grid of dynamical models to the ob- 
servations of Kepler, we have generated synthetic X-ray 
spectra for the emission from shocked ejecta in our hydro- 
dynamic simulations follow ing the methods presented in 
(Bad enes et al. | [200l |2005), using updated atomic data 
( Badenes et al.l 120061 ) . and in cluding radiati v e and ion- 
ization losses as described in iBadenes et al.l (|2007f) . In 
a constant density ambient medium, the synthetic spec- 
tra are controlled by three variables: the density of the 
ambient material, the SNR age, and the amount of col- 
lisionless heating between electrons and ions (/3) at the 
reverse shock. For our models where p a mb oc r~ 2 , the 
X-ray emission is instead influenced by the wind normal- 
ization (which sets the density). Higher values of M / (An 
fwind) will result in more circumstellar material closer 
to the explosion. We expect that this will lead to more 
collisional ionization in the shocked ejecta which will be 
directly reflected in the X-ray spectrum. 

The ability of our synthetic spectra to reproduce the 
observations is limited by the quality of the atomic data 
used in the code. The data for Ka blends are reasonably 
complete however, and we thus focus on those lines. As 
shown in Table [I] we have measured the line centroids 
and fluxes of several Ka lines. These can be compared 
directly against the synthetic spectra from our models. 
We focus only on the brightest lines of silicon, sulfur, and 
iron, and note the following additional constraints on our 
models- (1): the presence of Fe K emission at ~ 6.4 keV; 



(2): the absence of significant Si XIV Lya emission at 
~ 2.0 keV; (3): the amount of Fe L emission below 1 
keV, relative to oxygen suggests that the emission there 
is in fact ejecta and n ot shocked circumstellar material 
(|Revnolds et al.ll2007| ). 

In Figure [4] (left, upper panel) we plot the X-ray spec- 
trum from Figure [5] (left) compared against a subset of 
our spectral models. We plot the spectra from both the 
DDTa and DDTg models for a range of mass-loss rates, 
assuming wind velocities of 20 km s . Given constraints 
(1) and (3) above, we can rule out the DDTg models since 
they do not produce enough iron to match the emission 
seen both at 6.4 keV and Fe L emission at ~ 1 keV. By 
ruling out the DDTg models, we are also essentially rul- 
ing out the distance estimate of 3 < Dsnr < 6.4 kpc, 
if the blastwave is expanding into a p am b cx r~ 2 wind. 
In contrast, some DDTa models do reproduce the Fe Ka 
and Fe L emission. Those DDTa models with with M < 
4xl0 -6 M Q yr -1 are ruled out since they produce very 
little emission at 6.4 keV, in conflict with the data. 

In Figured] (left, upper panels), we plot the line cen- 
troids from the DDTa and DDTg models for various 
mass-loss rates for wind speeds of 20 km s _1 . The 
hatched region in each panel corresponds to the 90% con- 
fidence on the line centroid from the spectral fitting in 
Table [TJ While we have previously ruled out the DDTg 
models in this wind scenario based on the absence of Fe 
emission, we are also forced to rule out those DDTa mod- 
els that do produce Fe K emission as the predicted line 
centroids for nearly all models considered are at higher 
energies than what is measured. 

Finally, and for comparison with other possible CSM 
models, we note that those DDTa models shown in Fig- 
ure [4] that produce Fe K emission also result in signifi- 
cant Si XIV Lya emission at ~ 2.0 keV. This emission is 
not observed in the spectrum (Fig. [2]), and our spectral 
fits only resulted in an upper limit on the Si XIV Lya 
flux (Tabic Q}. This is shown quantitatively in Figured] 
(right) where we plot the ratio of the modeled flux from 
Si Ka to that from Si XIV Lya, for comparison against 
the measured ratio. In summary, based on the results of 
our comparisons between the observed and modeled bulk 
properties of the X-ray spectrum from the shocked ejecta, 
we must conclude that the X-ray emission from Kepler's 
SNR is not consistent with what is expected from su- 
pernova ejec ta interacting with a p a mb cx r~ 2 wind as 
proposed by iChiotellis etaLl (|2012D . 

3.3. Spectral Fitting: Wind + Cavity Model 

As discussed above, pure wind models do not provide 
a good fit to the X-ray spectrum. The reason is that 
the low wind velocities and moderate mass-loss rates 
that are required to match the SNR dynamics result in 
a significant amount of circumstellar material close to 
the explosion. Physically, this results in a stronger re- 
verse shock being driven into the expanding ejecta, re- 
sulting in higher average charge states- this affects both 
the Ka line centroid as well as the amount of Lya emis- 
sion. These are two measurable quantities which do not 
agree with the simple models described above. 

The possibility exists that the circumstellar envi- 
ronment is more complicated th an that of a simple 
r wind. Recurrent novae (e.g., iSokoloski et "all [2006; 
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iWood-Vasev fc Sokolo ski 2006) o r opt ically thick accre- 
tion winds (e.g.. lHachisu et al.l I1996D from the white 
dwarf could carve out a cavity around the progenitor, 
even while the companion is losing mass to a slow and 
dense wind. We consider the formation of a small (~ 
10 17 cm) cavity inside the dense wind. The cavity could 
form in a number of ways (see iChomiuk et~aLl 120121 
for a recent review), and time- variable Na I D absorp- 
tion spectra indicate that cavities can have sizes > 10 17 
cm (IPatat et al.l 120071 IBloirdin et al.l 120091: ISimon et all 
I2009t ISternberg et al.H201ll) . A small cavity such as this 
will not significantly alter the dynamics of the forward 
shock, but it will slow the ionization of the ejecta, re- 
sulting in lower charge states and dynamically younger 
shocked ejecta. 

In the bottom panel of Figure 0] (left), we plot the 
synthesized X-ray emission from models where the wind 
has a small cavity carved out of it. Additionally, we show 
in the bottom panel of Figure [5] (left) the computed line 
centroids for these models. In the right hand panel of 
Figure [5j we show the modeled flux ratio of Si Ka to 
Si XIV Lya compared against the upper limit on the flux 
ratio from the X-ray spectrum, as a function of mass- 
loss rate. As shown in these plots, the line centroids (Si 
Ka and S Ka) now agree with the measured values (for 
velocities < 20 km s _1 ). Additionally, these modified 
models show little Si XIV Lya emission at ~ 2 keV, as 
evidenced by the low flux ratios seen in the right hand 
panel, particularly for mass- loss rates less than 6xl0 -6 
M yr -1 . We can thus conclude that these models are 
consistent with the observed X-ray spectrum. 

When working within the dynamical constraints im- 
posed by the DDTa models, we require a distance greater 
than 7 kpc, and a wind luminosity greater than ~ 6x 10 32 
erg s" 1 . Additionally, we require iron emission, agree- 
ment between the modeled and measured Ka line cen- 
troids, and little or no Si XIV Lya emission. The high 
mechanical wind luminosity rules out models with low 
values of M, while the lower limit on the distance rules 
out the high mass-loss rate models. The presence of 
Si XIV Lya emission in the modeled spectra rule out 
those models with M > 6xl0~ 6 M Q yr _1 . We are left 
to conclude that a wind with a mass-loss rate of ~ 4- 
6 xl0~ 6 M Q yr~ 4 and wind velocities ~ 20 km s~ 4 
sculpted the CSM, and that a small cavity likely sur- 
rounded the progenitor, possibly caused by either recur- 
rent novae or optically thick accretion winds that pre- 
ceded the supernova. 

3.4. Spectral Fitting: Constant Density Models 

Finally, we consider models with a uniform a mbient 
mediu m, with density p a mb = 1-5 x 10 -24 g cm -3 . iVinkl 
(2008) measured the expansion of Kepler's SNR in the 
south and compared it against models for the expan- 
sion of a SNR blastwave in a constant density exter- 
nal medium. They found that models where p c) oc v~ 7 
were able to reproduce the kinematics. However, they 
concluded that the kinematics of the blastwave in the 
south point to an explosion with Esn < 5 x 10 50 erg 
for the canonical distance of 4 kpc, assuming an ambi- 
ent medium density of nu = 1 cm~ 3 . Higher explosion 
energies work if the SNR is at a distance considerably 
greater than 4 kpc, or if the ambient medium density > 



10 cm' 3 . 

In Figure \$\(right), we show the computed X-ray emis- 
sion for DDTa models in a constant density medium. As 
shown in Fig. 2J these models qualitatively reproduce 
the bulk features seen in the spectrum, and an ambi- 
ent medium density of 2-5 x 10~ 24 g cm~ 3 seems to be 
able to reproduce both the line centroids as well as the 
line flux ratios. This is confirmed in Figure [6j where 
we see that models with densities of < 3xl0~ 24 g cm -3 
produce the required Si XIV Lya emission and also pre- 
dict line centroids in agreement with the measured val- 
ues. For these models, the forward shock radii at 400 
yrs range from 2.9 - 3.8 pc. For the measured angu- 
lar radius of 2', this corresponds to a distance to Kepler 
of 5 - 6.5 kpc. cons istent with the distance estimate of 
ISankrit et all ()2005h . 

While the constant density models can explain the ob- 
served X-ray emission and also provide estimates for the 
distance to Kepler that are in agreemen t with proper 
motio n measurements of Balmer filaments (|Sankrit et al.l 
2005), it is difficult to reconcile a constant density CSM 
with the bowshock structure in the north. To form a 
bowshock, a converging flow is required between an out- 
flow from the progenitor and the bulk motion of the am- 
bient ISM (in the rest frame of the progenitor). Two 
alternative scenarios are considered: First, for either a 
wind that is not spherically sym metric but consists o f 
equatorial and polar components (jOwocki et al.l [19941) . 
or other non-sp herical outflows , such as in the recurrent 
nova RS Oph (ISokolosk i et al.|[200l IDrake et aL| [200l 
lOrlando et al.ll2009D . the CSM could be sufficiently mod- 
ified such that it is better described as a constant density 
environment, rather than a r -2 wind. 

Alternatively, we can consider that Kepler's SNR 
exploded in an environment with a density gradient. 
iPetrukl (|1999j ) modeled the evolution of a SNR in a large 
density gradient. In Figure 5(a) of IPetrukl (fl999h . the 
surface brightness profiles for the modeled remnant are 
qualitatively similar to the morphology of Kepler. That 
is, a bowshock like structure forms i n the d i rectio n of 
the density gradient. For RCW 86, IPetrukl (|1999| ) re- 
quire a density contrast of < 2 between the north and 
south, in order to reproduce the observed X-ray emis- 
sion at an age of 1800 yr. For Kepler, a density con- 
trast considerably higher than 2 would be required to 
produce the same morpholo gy in just 400 yr. 160/xm 
emission (jGomez et al.l l2012) suggests a density gradient 
towards the north-northwest of Kepler, but they point 
out that much of that e mission may arise f rom foreground 
or background sources. (jBlair et al.ll2007|) reached a sim- 
ilar conclusion from Spitzer 160/im observations, noting 
a north-south gradient in the emission at that wave- 
length, indicative of a density gradient in the circum- 
stellar medium. If the bowshock structure is caused by 
a gradient in the density of the circumstellar environ- 
ment, then that would suggest that the SN exploded in 
an environment that was closer to uniform in density, 
rather than one shaped by a slow and dense stellar wind 
interacting with the ISM. If this were indeed the case 
for Kepler, then it would be similar to the candidate 
Type la remnant G299.2-2.9. In the case of that rem- 
nant, a signficant density gradient (3x) is required to 
expl ain the observed asymmetry in the X-ray morphol- 
ogy ()Park et al.ll2007D . However, that SNR is estimated 
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to be ten times older than Kepler, so the spatial scale of 
the density gradient is correspondingly larger. 

4. CONCLUSIONS 

In order to explain the observed dynamics and X-ray 
emission of Kepler's supernova remnant, we have coupled 
models for ejecta from Type la SNe to a model for a 
circumstellar medium shaped by a slow stellar wind. We 
have synthesized the X-ray emission from these models 
and compared it to the observed X-ray spectrum. Our 
conclusions are: 

1. If the distance to Kepler is taken to be less than 
6.4 kpc, then the dynamics (e.g. the forward shock 
radius) are better explained by a subenergetic ex- 
plosion model (DDTg) interacting with a slow and 
dense stellar wind. However, while the subener- 
getic models can explain the dynamics, the syn- 
thesized X-ray spectra do not agree with the ob- 
served X-ray emission. Specifically, the subener- 
getic model does not produce enough iron to ex- 
plain the observed emission. The presence of strong 
Fe emission in the spectrum rules out all DDTg 
models. 

2. If the distance is taken to be greater than 7 kpc, 
then the dynamics can be explained by the ener- 
getic la models. While the X-ray spectrum of Ke- 
pler shows significant emission from Fe K and Fe 
L, there is little emission from H-like Lya of inter- 
mediate mass elements such as silicon. The mod- 
eled line centroids in the DDTa + wind models do 
not agree with the measured line centroids, even 
in those models with high mechanical luminosities 
(i.e. high wind speeds). Additionally, the pure 
wind models over predict the amount of Si XIV 
Lya emission relative to that which is observed. 
This rules out the DDTa + wind models with CSM 
winds that extend all the way to the explosion site. 

3. Since the DDTg models are ruled out based on the 
presence of Fe emission, and the DDTa + wind 
models are ruled out based on the measured K- 
shell line centroids and Si Ka to Si XIV Lya flux 
ratio, we consider the presence of a small, low den- 
sity cavity around the progenitor prior to the SN. 
The low density cavity alters the ionization state of 
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the shocked ejecta but does not significantly change 
the dynamics of the forward shock. We find that 
the presence of the cavity reduces both the amount 
of Si XIV Lya emission as well as the modeled line 
centroids to values consistent with the measured 
centroids. We find that the energetic la model 
(DDTa) interacting with a wind with velocity < 20 
km s _1 and mass-loss rate of < 6xl0 -6 M Q yr" 1 
is a good match both spectrally and dynamically, 
assuming a distance > 7 kpc. 

4. DDTa models with a constant density ambient 
medium (p am b « 2 x 10~ 24 g cm~ 3 ) can explain 
the observed X-ray emission and place the SNR 
at a dist ance that is consi stent with previous 
estimates (jSankrit et al.ll2005T) . In order to explain 
the bowshock structure in the north, a north-south 
density gradient in the interstellar medium is 
required. Evidence for such a gradient is seen in 
160/im emission maps of Kepler and its vicinity 
([Blair et al.M2007h . 

Thus, we conclude that in order to explain both the 
dynamics and X-ray emission of Kepler in the context 
of a wind environment, the explosion was energetic, pro- 
ducing ~ 1M Q of 56 Ni. The supernova occurred in a 
slow and dense CSM wind with a small central cavity. 
The observed X-ray emission, which traces the explo- 
sion energetics, requires a distance > 7 kpc, in the CSM 
wind model. Alternatively, we find that the X-ray emis- 
sion and dynamics of Kepler can also be explained by 
a constant CSM density model. These models place the 
SNR at a distance of ~ 5-6 kpc, consistent with previous 
measurements. In the constant density scenario, a north- 
south density gradient in the ISM would be required in 
order to explain the SNR morphology. 
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Fig. 1. — Kepler's SNR viewed in X-rays with Chandra ACIS-S3. The RGB image shows 0.4-0.75 keV emission in red, 0.75-1.2 keV 
emission in green, and 1.2—7.0 keV emission in blue. In the image, north is up and east is to the left. 
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Fig. 2. — Left: Chandra ACIS-S spectrum of the pie-sliced region of Kepler's SNR shown in Figure [T] For reference, we mark the 
positions of several Ka lines that are used in the analysis. Right: Spectral fit to the 1.7-7.0 keV spectrum of the pie-sliced region shown in 
Figure[T]The fit parameters are listed in Table[T] and the positions of the Ka lines are also marked. In both figures, the data are unbinned 
(1024 PHA channels). 
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Fig. 3. — Forward shock radius as a function of wind luminosity for a range of mass-loss rates and wind speeds, for Type la SN ejecta 
models DDTa and DDTg. On the right hand axis is the distance to Kepler that is implied by the forward shock radius, for an angular 
diameter of 4', but the models may be scaled to any distance. 
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Fig. 4. — Left: Comparisons between the modeled X-ray spectrum and the spectrum shown in Figure[T] In the upper panel, we compare 
the DDTa (solid lines) and DDTg (dashed lines) modeled spectra for a range of mass-loss rates, assuming a wind speed of 20 km s~ 1 . In 
the bottom panel, we plot the DDTa models for a wind speed of 20 km s _1 and a range of mass-loss rates, assuming a small, low density 
cavity has formed around the progenitor. Bight: DDTa models in a constant density ambient medium. 



Wind Models 



5. 0.8 
? 0.6 

6 

— 0.4 
■3 

0.2 




I 1 S 1 I ' 
= S Ka 



I 1^ I I I I I I I I I j I I I I L I I 

= FeKa 



* 



1.84 1.845 1.85 2.35 2.4 2.45 6.45 6.5 6.55 6.6 6.65 

Centroid [keV] 



Wind + Cavity Models 




10 km s-' = S Ka 
• 20 km s- 1 



■ . •■ moss- 
loss rates 



I I I i I I I I i I I 




1.84 1.845 1.85 2.35 2.4 2.45 6.45 6.5 6.55 6.6 6.65 

Centroid [keV] 



DDTo (ft « r-2) 
20 km s"' 

fl * DDTg (p « r" 2 ) 

„ = 20 km s"' 
DDTo (Cavity) 
„ = 10 km s"' 
DDTa (Cavity) 
= 20 km s-' 




favored by 
data 



1 1,5 2 

log,„[Si Ka/Si Lya] 



Fig. 5. — Left: In the upper panel, we plot the measured versus computed line centroids for Si Ka, S Ka, and Fe Ka, for the DDTa 
(black solid) and DDTg (red dashed) models. The hatched region in each panel corresponds to the measured centroid including the 90% 
confidence interval The DDTg models and a subset of the DDTa models (M < 4xl0 — 6 Mq yr _1 ) do not produce any Fe K emission. In 
the lower panel, we plot the line centroids for the DDTa cavity models for v w ; nt j = 10 (green solid) and 20 (blue dashed) km s —1 , for a 
range of mass-loss rates. The line centroids and errors are indicated by the vertical hatched regions in each panel. The allowed mass-loss 
rates as dictated by the comparison between the measured and modeled line centroids are marked by the horizontal cross hatched region. 
Right: Si Ka/Si XIV Lya flux ratio from each model as compared against the measured value from Table [T] The flux ratios are plotted 
as a function mass-loss rate for both the pure wind models as well as the models that include a cavity. The line styles are the same as in 
the left hand plot. We mark the region of the plot that is not consistent with the measured line ratio which rules out almost all pure wind 
models. 
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Fig. 6. — Top: Si Ka/Si XIV Lya flux ratio for energetic and subenergetic models in a constant density environment, for either 
instantaneous (J3 = 1) or mass proportional (/3 = 0.1) electron heating. The modeled ratio is relatively insensitive to the electron heating 
model. Bottom: Measured versus computed line centroids for the DDTa (blue) and DDTg (red) models for Si Ka, S Ka, and Fe Ka. 
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TABLE 1 

Spectral Parameters in Kepler's SNR 





Parameter 


Fitted Valuc a 




Power-Law Continuum 


a 

Norm (lCT 3 


photons cm 2 s 1 at 1 keV) 


9 R7 +0.06 
Z - D < -0.01 
1 90+ ' 20 


Line Fluxes (10 ° photons cm 




Si Ka b 
Si XIV Lya 
Si K/3 
S Ka 
Ar Ka 
Ca Ka 
Fe Ka 




49.9±0.20 
< 4.70±0.47 
5.64±0.30 
13.5±0.30 
1.19±0.16 
0.38±0.07 
3.84±0.14 


Line Centroids (keV) 


Si Ka 
Si K/3 
S Ka 
Ar Ka 
Ca Ka 
Fe Ka 




9 1Qn+ 004 

2.iyu_ 003 

9 /I9C+0.007 
°— 0.006 

o n77 +o.oos 

°- u ' ' -0.006 
o 7QQ+0.010 
3- '99-0.010 

fi 45n +0010 


Line Width (eV) 


Si Ka 
Si K/3 
S Ka 
Ar Ka 
Ca Ka 
Fe Ka 




25.9l^ 
35.91^ 

— 5.5 

83.01*;* 



a Errors correspond to 90% confidence intervals 

b Ka emission here refers to all allowed transitions to the 

ground state 



